*****************************************************************
*****************************************************************
*****                                                       *****
*****       Mona Morgan-Collins (Durham University, UK)     *****
*****       Grace Natusch (The Civil Service, UK)           *****
*****       Contact: mona.morgan-collins@durham.ac.uk       *****
*****                                                       *****
*****       At the Intersection of Gender and Class:        *****
***** How Were Newly Enfranchised Women Mobilized in Sweden?*****
*****                                                       *****
*****           Comparative Political Sudies                *****
*****                                                       *****
*****         Replicating Analyses in the Appendix          *****
*****                                                       *****
*****************************************************************
*****************************************************************


*****************************************************
*Figure A1: Södertälje's Neighborhoods and Housing
*****************************************************

*Nothing to replicate


*****************************************************
*Figure A2: Class by Neighborhood and Year 
*****************************************************
use dta\sod2134, clear                                           //using sod2134 data set

*Sub-Figure c: Municipal 1921
#delimit ; 
graph bar worker upmid if year==1921, over(neighb  , relabel(1 "Al" 2 "At"  3 "Br" 4 "La" 5 "Ro" 6 "St R" 7"T") 
label(labsize(huge))) blabel(bar, size(large) format(%9.2f)) ytitle("Class (%)", size(huge))  ylabel (0(0.25)1 , labsize(huge) nogrid) 
box(1, color(gray)) box(2, color(gs4)) marker(1 , mcolor(gray)) title("",  size(vlarge))  legend(off) scheme(s1color) ysize(5) ; 
#delimit cr

*Sub-Figure d: Municipal 1934
#delimit ; 
graph bar worker upmid if year==1934, over(neighb  , relabel(1 "Al" 2 "At"  3 "Br" 4 "La" 5 "Ro" 6 "St R" 7"T") 
label(labsize(huge))) blabel(bar, size(large) format(%9.2f)) ytitle("Class (%)", size(huge))  ylabel (0(0.25)1 , labsize(huge) nogrid) 
box(1, color(gray)) box(2, color(gs4)) marker(1 , mcolor(gray)) title("",  size(vlarge))  legend(off) scheme(s1color) ysize(5) ; 
#delimit cr


*****************************************************
*Table A1: Comparing Södertälje with Sweden
*****************************************************

*Nothing to replicate


*****************************************************
*Table A2: Summary Statistics of Eligibles
*****************************************************
use dta\sod2134, clear                                           //using sod2134 data set

sum voted_munic if year==1921 & right_munic==1                   //turnout by elecions
sum voted_county if year==1921 & right_county==1
sum voted_munic if year==1934 & right_munic==1

sum female if year==1921 & right_munic==1                        //sex by election
sum female if year==1921 & right_county==1
sum female if year==1934 & right_munic==1

sum age if year==1921 & right_munic==1                           //age by election
sum age if year==1921 & right_county==1
sum age if year==1934 & right_munic==1

sum married if year==1921 & right_munic==1                       //married by election
sum married if year==1921 & right_county==1
sum married if year==1934 & right_munic==1

sum worker if year==1921 & right_munic==1                        //working class by election
sum worker if year==1921 & right_county==1
sum worker if year==1934 & right_munic==1

sum upmid if year==1921 & right_munic==1                         //upper and middle class by election
sum upmid if year==1921 & right_county==1
sum upmid if year==1934 & right_munic==1


*****************************************************
*Table A3-A5: Classification of Occupations (Part I-3)
*****************************************************
use dta\sod2134, clear                       //using sod2134 data set

tab occup_eng if occup_cat=="I"      //agriculture - unclassified
tab occup_eng if occup_cat=="I_a"    //agriculture - livestock & farm
tab occup_eng if occup_cat=="I_b"    //agriculture - fishery
tab occup_eng if occup_cat=="I_c"    //agriculture - forestry
tab occup_eng if occup_cat=="II"     //industry - unclassified
tab occup_eng if occup_cat=="II_a"   //industry - food
tab occup_eng if occup_cat=="II_b"   //industry - textile
tab occup_eng if occup_cat=="II_c"   //industry - wood
tab occup_eng if occup_cat=="II_d"   //industry - construction
tab occup_eng if occup_cat=="II_e"   //industry - leather, hair & rubber
tab occup_eng if occup_cat=="II_f"   //industry - stone
tab occup_eng if occup_cat=="II_g"   //industry - mining & metal
tab occup_eng if occup_cat=="II_h"   //industry - paper & graphics
tab occup_eng if occup_cat=="II_i"   //industry - chemical
tab occup_eng if occup_cat=="III"    //transport & commmunication
tab occup_eng if occup_cat=="IV"     //unclassified Manual
tab occup_eng if occup_cat=="V"      //domestic service
tab occup_eng if occup_cat=="VI_a"   //public sector - administration, defence
tab occup_eng if occup_cat=="VI_b"   //public sector - religion, teaching, art
tab occup_eng if occup_cat=="VI_c"   //public sector - poor, health, social worker
tab occup_eng if occup_cat=="VII"    //sales and clerical sector
tab occup_eng if occup_cat=="VIII"   //merchants, owners and tenants
tab occup_eng if occup_cat=="IX"     //professions
tab occup_eng if occup_cat=="X"      //large owners
tab occup_eng if occup_cat=="0"      //students, dependents


******************************************************************
*Figure A3: Turnout by Gender, Class & Election in Södertälje
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Municipal 1921
#delimit ;
spineplot votedm_rec group if year==1921, xlabel(, angle(v) axis(2) labsize(vlarge)) xlabel(, axis(1) labsize(vlarge))  ylabel(, axis(2) labsize(vlarge)) 
xtitle("", axis(2) size(vlarge)) xtitle("Fraction by Class and Sex", axis(1) size(vlarge)) ytitle("Fraction by Voted", axis(2) size(vlarge)) ytitle("", axis(1) size(vlarge))
bar1(bcolor(gs15)) bar2(bcolor(gray)) scheme(s1mono) legend(off) ysize(6); 
#delimit cr

*Sub-Figure a: County 1921
#delimit ;
spineplot votedc_rec group if year==1921, xlabel(, angle(v) axis(2) labsize(vlarge)) xlabel(, axis(1) labsize(vlarge))  ylabel(, axis(2) labsize(vlarge)) 
xtitle("", axis(2) size(vlarge)) xtitle("Fraction by Class and Sex", axis(1) size(vlarge)) ytitle("Fraction by Voted", axis(2) size(vlarge)) ytitle("", axis(1) size(vlarge))
bar1(bcolor(gs15)) bar2(bcolor(gray)) scheme(s1mono) legend(off) ysize(6); 
#delimit cr

*Sub-Figure c: Municipal 1934
#delimit ;
spineplot votedm_rec group if year==1934, xlabel(, angle(v) axis(2) labsize(vlarge)) xlabel(, axis(1) labsize(vlarge))  ylabel(, axis(2) labsize(vlarge)) 
xtitle("", axis(2) size(vlarge)) xtitle("Fraction by Class and Sex", axis(1) size(vlarge)) ytitle("Fraction by Voted", axis(2) size(vlarge)) ytitle("", axis(1) size(vlarge))
bar1(bcolor(gs15)) bar2(bcolor(gray)) scheme(s1mono) legend(off) ysize(6); 
#delimit cr
 

******************************************************************
*Figure A4: Plotting Class Composition of neighbors on Property Size
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Workers 1921
#delimit ; 
 binscatter tworkp_pc propsize if year==1921, lc(black) mc(gray) m(Oh) n(10) 
 ytitle(Workers in a Property (%), size(vlarge)) ylabel(0(25)100, labsize(vlarge)) xtitle(# Electors per Property (Binned by Decile), size(large)) xlabel(0(25)100, labsize(vlarge)) 
 legend(off) scheme(s1color) ysize(6);
#delimit cr

*Sub-Figure b: Upper&Middle Class 1921
#delimit ; 
 binscatter tupmidp_pc propsize if year==1921, lc(black) mc(gray) m(Oh) n(10) 
 ytitle(Upper & Middle Class in a Property (%), size(vlarge)) ylabel(0(25)100, labsize(vlarge)) xtitle(# Electors per Property (Binned by Decile), size(large)) xlabel(0(25)100, labsize(vlarge)) 
 legend(off) scheme(s1color) ysize(6);
#delimit cr

*Sub-Figure c: Workers 1934
#delimit ; 
 binscatter tworkp_pc propsize if year==1934, lc(black) mc(gray) m(Oh) n(10) 
 ytitle(Workers in a Property (%), size(vlarge)) ylabel(0(25)100, labsize(vlarge)) xtitle(# Electors per Property (Binned by Decile), size(large)) xlabel(0(25)100, labsize(vlarge)) 
 legend(off) scheme(s1color) ysize(6);
#delimit cr

*Sub-Figure d: Upper&Middle Class 1934
#delimit ; 
 binscatter tupmidp_pc propsize if year==1934, lc(black) mc(gray) m(Oh) n(10) 
 ytitle(Upper & Middle Class in a Property (%), size(vlarge)) ylabel(0(25)100, labsize(vlarge)) xtitle(# Electors per Property (Binned by Decile), size(large)) xlabel(0(25)100, labsize(vlarge)) 
 legend(off) scheme(s1color) ysize(6);
#delimit cr
 

*****************************************************
*Table A6: Turnout by Gender, Class and Election
*****************************************************
use dta\sod2134, clear                       //using sod2134 data set

sum voted_munic if year==1921 & right_munic==1                  //turnout by election
sum voted_county if year==1921 & right_county==1
sum voted_munic if year==1934 & right_munic==1

sum voted_munic if year==1921 & right_munic==1 & female==1      //turnout by election and sex
sum voted_county if year==1921 & right_county==1 & female==1
sum voted_munic if year==1934 & right_munic==1 & female ==1
sum voted_munic if year==1921 & right_munic==1 & female==0
sum voted_county if year==1921 & right_county==1 & female==0
sum voted_munic if year==1934 & right_munic==1 & female ==0

sum voted_munic if year==1921 & right_munic==1 & worker==1      //turnout by election and class
sum voted_county if year==1921 & right_county==1 & worker==1
sum voted_munic if year==1934 & right_munic==1 & worker==1
sum voted_munic if year==1921 & right_munic==1 & upmid==1
sum voted_county if year==1921 & right_county==1 & upmid==1
sum voted_munic if year==1934 & right_munic==1 & upmid==1


*****************************************************
*Table A7: Observed Demographics by Class and Election
*****************************************************
use dta\sod2134, clear                       //using sod2134 data set

sum married  female age  if right_munic==1   & year==1921  & worker==1  
sum married  female age  if right_munic==1   & year==1921  & upmid==1  
sum married  female age  if right_county==1  & year==1921  & worker==1
sum married  female age  if right_county==1  & year==1921  & upmid==1
sum married  female age  if right_munic==1   & year==1934  & worker==1
sum married  female age  if right_munic==1   & year==1934  & upmid==1


*****************************************************
*Table A8: Comparing Single and Multi-Family Properties
*****************************************************
use dta\sod2134, clear                       //using sod2134 data set
                                       
sum married female age worker upmid   if sinfampro==1 & year==1921               
sum married female age worker upmid   if sinfampro==0 & year==1921
sum married female age worker upmid   if sinfampro==1 & year==1934
sum married female age worker upmid   if sinfampro==0 & year==1934


******************************************************************
*Figure A5: Average Neighbor Effects for Class and Gender
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*sub-figure a: proportion of women workers by gender and class composition of properties
#delimit ; 
twoway hist tworkfp_pc if year==1921 & class_wife==3 & female==1, lcolor(black) color(black) percent w(2)  || hist tworkfp_pc if year==1934 & class_wife==3 & female==1, lcolor(gray) color(none) percent w(2)  ||
, ytitle(% Working class women, size(vlarge)) ylabel(0(10)30 , labsize(vlarge))  xtitle("% Working class women neighbors", size(vlarge)) xlabel(0(20)100, labsize(vlarge))  
legend(off) scheme(s1color) ysize(5) ; 
#delimit cr

*Sub-Figure b: Average Neighbor Effects by Gender and Class Composition of Properties 
gen f21mw = workXtworkfp_pc
gen f21cw = workXtworkfp_pc
gen f34mw = workXtworkfp_pc

quietly reg voted_munic  i.worker f21mw   age age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
quietly estimates store f21mw
quietly reg voted_county i.worker f21cw   age age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
quietly estimates store f21cw
quietly reg voted_munic  i.worker f34mw   age age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
quietly estimates store f34mw

#delimit ;                                      
coefplot 
(f21mw, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(f21cw, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(f34mw, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)4)) 
keep(  f21mw f21cw f34mw  )  coeflabels(  f21mw="M21" f21cw="C21" f34mw="M34"        )  
yline(0, lcolor(gray))  levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5)  ;
#delimit cr


******************************************************************
*Table A9: Average Neighbor Effects By Class and Gender  (Full Results for Working-Class)
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

 reg voted_munic  i.worker w21m    C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
 reg voted_county i.worker w21c    C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
 reg voted_munic  i.worker w34m    C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
 reg voted_munic  i.worker w21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
 reg voted_county i.worker w21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
 reg voted_munic  i.worker w34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
 reg voted_munic  i.worker w21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
 reg voted_county i.worker w21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
 reg voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 
 

******************************************************************
*Table A10: Average Neighbor Effects By Class and Gender  (Full Results for Upper and Middle Class)
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

 reg voted_munic  i.upmid um21m    C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
 reg voted_county i.upmid um21c    C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
 reg voted_munic  i.upmid um34m    C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
 reg voted_munic  i.upmid um21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
 reg voted_county i.upmid um21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
 reg voted_munic  i.upmid um34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
 reg voted_munic  i.upmid um21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
 reg voted_county i.upmid um21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
 reg voted_munic  i.upmid um34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 


******************************************************************
*Figure A6: Alternative Measurement of Class Composition
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Working Class
gen w21m_r = workXtworkp_pc_r
gen w21c_r = workXtworkp_pc_r
gen w34m_r = workXtworkp_pc_r
gen w21mf_r = workXtworkp_pc_r
gen w21cf_r = workXtworkp_pc_r
gen w34mf_r = workXtworkp_pc_r
gen w21mm_r = workXtworkp_pc_r
gen w21cm_r = workXtworkp_pc_r
gen w34mm_r = workXtworkp_pc_r

quietly reg voted_munic  i.worker w21m_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
quietly estimates store w21m_r
quietly reg voted_county i.worker w21c_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
quietly estimates store w21c_r
quietly reg voted_munic  i.worker w34m_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
quietly estimates store w34m_r
quietly reg voted_munic  i.worker w21mf_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
quietly estimates store w21mf_r
quietly reg voted_county i.worker w21cf_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
quietly estimates store w21cf_r
quietly reg voted_munic  i.worker w34mf_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
quietly estimates store w34mf_r
quietly reg voted_munic  i.worker w21mm_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
quietly estimates store w21mm_r
quietly reg voted_county i.worker w21cm_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
quietly estimates store w21cm_r
quietly reg voted_munic  i.worker w34mm_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 
quietly estimates store w34mm_r
 
#delimit ;                                      
coefplot 
(w21m_r, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21c_r, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w34m_r, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21mf_r, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cf_r, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mf_r, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mm_r, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w21cm_r, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w34mm_r, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( w21m_r w21c_r w34m_r  w21mf_r w21cf_r w34mf_r  w21mm_r w21cm_r w34mm_r )  coeflabels(w21m_r="M21" w21c_r="C21" w34m_r="M34"   w21mf_r="M21" w21cf_r="C21" w34mf_r="M34"  w21mm_r="M21" w21cm_r="C21" w34mm_r="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 8 "Men", color(black) size(large)) text(-0.0045 5 "Women", color(gray) size(large)) text(-0.0045 2 "All", color(red) size(large)) ;
#delimit cr

*Sub-Figure b: Upper and Middle Class
gen um21m_r = upmidXtupmidp_pc_r
gen um21c_r = upmidXtupmidp_pc_r
gen um34m_r = upmidXtupmidp_pc_r
gen um21mf_r = upmidXtupmidp_pc_r
gen um21cf_r = upmidXtupmidp_pc_r
gen um34mf_r = upmidXtupmidp_pc_r
gen um21mm_r = upmidXtupmidp_pc_r
gen um21cm_r = upmidXtupmidp_pc_r
gen um34mm_r = upmidXtupmidp_pc_r

quietly reg voted_munic  i.upmid um21m_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
quietly estimates store um21m_r
quietly reg voted_county i.upmid um21c_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
quietly estimates store um21c_r
quietly reg voted_munic  i.upmid um34m_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
quietly estimates store um34m_r
quietly reg voted_munic  i.upmid um21mf_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
quietly estimates store um21mf_r
quietly reg voted_county i.upmid um21cf_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
quietly estimates store um21cf_r
quietly reg voted_munic  i.upmid um34mf_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
quietly estimates store um34mf_r
quietly reg voted_munic  i.upmid um21mm_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
quietly estimates store um21mm_r
quietly reg voted_county i.upmid um21cm_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
quietly estimates store um21cm_r
quietly reg voted_munic  i.upmid um34mm_r   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 
quietly estimates store um34mm_r
 
#delimit ;                                      
coefplot 
(um21m_r, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um21c_r, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um34m_r, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um21mf_r, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um21cf_r, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um34mf_r, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um21mm_r, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(um21cm_r, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(um34mm_r, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( um21m_r um21c_r um34m_r  um21mf_r um21cf_r um34mf_r  um21mm_r um21cm_r um34mm_r )  coeflabels(um21m_r="M21" um21c_r="C21" um34m_r="M34"   um21mf_r="M21" um21cf_r="C21" um34mf_r="M34"  um21mm_r="M21" um21cm_r="C21" um34mm_r="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 8 "Men", color(black) size(large)) text(-0.0045 5 "Women", color(gray) size(large)) text(-0.0045 2 "All", color(red) size(large)) ;
#delimit cr 


******************************************************************
*Figure A7: Alternative Measurement of Class Homogeneity
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

gen m21m = midXtmidp_pc
gen m21c = midXtmidp_pc
gen m34m = midXtmidp_pc
gen m21mf = midXtmidp_pc
gen m21cf = midXtmidp_pc
gen m34mf = midXtmidp_pc
gen m21mm = midXtmidp_pc
gen m21cm = midXtmidp_pc
gen m34mm = midXtmidp_pc

quietly reg voted_munic  i.mid m21m  C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
quietly estimates store m21m
quietly reg voted_county i.mid m21c   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
quietly estimates store m21c
quietly reg voted_munic  i.mid m34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
quietly estimates store m34m
quietly reg voted_munic  i.mid m21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
quietly estimates store m21mf
quietly reg voted_county i.mid m21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
quietly estimates store m21cf
quietly reg voted_munic  i.mid m34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
quietly estimates store m34mf 
quietly reg voted_munic  i.mid m21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
quietly estimates store m21mm
quietly reg voted_county i.mid m21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
quietly estimates store m21cm
quietly reg voted_munic  i.mid m34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 
quietly estimates store m34mm
 
#delimit ;                                      
coefplot 
(m21m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(m21c, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(m34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(m21mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(m21cf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(m34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(m21mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(m21cm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(m34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( m21m m21c m34m  m21mf m21cf m34mf  m21mm m21cm m34mm )  coeflabels(m21m="M21" m21c="C21" m34m="M34"   m21mf="M21" m21cf="C21" m34mf="M34"  m21mm="M21" m21cm="C21" m34mm="M34", angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 8 "Men", color(black) size(large)) text(-0.0045 5 "Women", color(gray) size(large)) text(-0.0045 2 "All", color(red) size(large)) ;
#delimit cr
 

******************************************************************
*Figure A8: Excluding Multiwife Properties
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Working Class
quietly reg voted_munic  i.worker w21m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921  &  wifemulti<=1
quietly estimates store w21m
quietly reg voted_county i.worker w21c   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921  &  wifemulti<=1
quietly estimates store w21c
quietly reg voted_munic  i.worker w34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  &  wifemulti<=1
quietly estimates store w34m
quietly reg voted_munic  i.worker w21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1  &  wifemulti<=1
quietly estimates store w21mf
quietly reg voted_county i.worker w21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1  &  wifemulti<=1
quietly estimates store w21cf
quietly reg voted_munic  i.worker w34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  &  wifemulti<=1
quietly estimates store w34mf 
quietly reg voted_munic  i.worker w21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  &  wifemulti<=1
quietly estimates store w21mm
quietly reg voted_county i.worker w21cm   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0  &  wifemulti<=1
quietly estimates store w21cm
quietly reg voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  &  wifemulti<=1
quietly estimates store w34mm
 
#delimit ;                                      
coefplot 
(w21m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21c, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w21cm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( w21m w21c w34m  w21mf w21cf w34mf  w21mm w21cm w34mm )  coeflabels(w21m="M21" w21c="C21" w34m="M34"   w21mf="M21" w21cf="C21" w34mf="M34"  w21mm="M21" w21cm="C21" w34mm="M34", angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 8 "Men", color(black) size(large)) text(-0.0045 5 "Women", color(gray) size(large)) text(-0.0045 2 "All", color(red) size(large)) ;
#delimit cr

*Sub-Figure b: Upper and Middle Class
quietly reg voted_munic  i.upmid um21m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921  &  wifemulti<=1  
quietly estimates store um21m
quietly reg voted_county i.upmid um21c   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921  &  wifemulti<=1
quietly estimates store um21c
quietly reg voted_munic  i.upmid um34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  &  wifemulti<=1
quietly estimates store um34m
quietly reg voted_munic  i.upmid um21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1  &  wifemulti<=1
quietly estimates store um21mf
quietly reg voted_county i.upmid um21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1  &  wifemulti<=1
quietly estimates store um21cf
quietly reg voted_munic  i.upmid um34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  &  wifemulti<=1
quietly estimates store um34mf 
quietly reg voted_munic  i.upmid um21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  &  wifemulti<=1
quietly estimates store um21mm
quietly reg voted_county i.upmid um21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0  &  wifemulti<=1
quietly estimates store um21cm
quietly reg voted_munic  i.upmid um34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  &  wifemulti<=1
quietly estimates store um34mm
 
#delimit ;                                      
coefplot 
(um21m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um21c, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um21mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um21cf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um21mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(um21cm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(um34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( um21m um21c um34m  um21mf um21cf um34mf  um21mm um21cm um34mm )  coeflabels(um21m="M21" um21c="C21" um34m="M34"   um21mf="M21" um21cf="C21" um34mf="M34"  um21mm="M21" um21cm="C21" um34mm="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 8 "Men", color(black) size(large)) text(-0.0045 5 "Women", color(gray) size(large)) text(-0.0045 2 "All", color(red) size(large)) ;
#delimit cr 


******************************************************************
*Figure A9: Alternative Specification - Excluding Large Properties in 1934
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Working Class
quietly reg voted_munic  i.worker w34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  & propsize<25
quietly estimates store w34m
quietly reg voted_munic  i.worker w34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & propsize<25
quietly estimates store w34mf 
quietly reg voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & propsize<25
quietly estimates store w34mm
 
#delimit ;                                      
coefplot 
(w34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(1(1)3)) 
keep(   w34m    w34mf    w34mm )  coeflabels( w34m="M34"    w34mf="M34"   w34mm="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(1.5 2.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 3 "Men", color(black) size(large)) text(-0.0045 2 "Women", color(gray) size(large)) text(-0.0045 1 "All", color(red) size(large)) ;
#delimit cr

*Sub-Figure b: Upper and Middle Class
quietly reg voted_munic  i.upmid um34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  & propsize<25
quietly estimates store um34m
quietly reg voted_munic  i.upmid um34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & propsize<25
quietly estimates store um34mf 
quietly reg voted_munic  i.upmid um34mm   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & propsize<25
quietly estimates store um34mm
 
#delimit ;                                      
coefplot 
(um34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(1(1)3)) 
keep(   um34m    um34mf    um34mm )  coeflabels( um34m="M34"    um34mf="M34"   um34mm="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(1.5 2.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 3 "Men", color(black) size(large)) text(-0.0045 2 "Women", color(gray) size(large)) text(-0.0045 1 "All", color(red) size(large)) ;
#delimit cr
 

******************************************************************
*Figure A10: Excluding One Neighborhood at the Time
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

foreach i in 1 2 3 4 5 6 7    {
 quietly reg voted_munic  i.worker w21m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921     & n`i'==0  
 quietly estimates store w21m`i'
 quietly reg voted_county  i.worker w21c   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1  & year==1921   & n`i'==0  
 quietly estimates store w21c`i'  
 quietly reg voted_munic  i.worker w34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934     & n`i'==0  
 quietly estimates store w34m`i'
 quietly reg voted_munic  i.worker w21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1  & n`i'==0  
 quietly estimates store w21mf`i'
 quietly reg voted_county  i.worker w21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1  & year==1921   & female==1  & n`i'==0  
 quietly estimates store w21cf`i'  
 quietly reg voted_munic  i.worker w34mf   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & n`i'==0  
 quietly estimates store w34mf`i'
 quietly reg voted_munic  i.worker w21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  & n`i'==0  
 quietly estimates store w21mm`i'
 quietly reg voted_county  i.worker w21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1  & year==1921   & female==0  & n`i'==0  
 quietly estimates store w21cm`i'  
 quietly reg voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & n`i'==0  
 quietly estimates store w34mm`i' 
 } 

foreach i in 1 2 3 4 5 6 7    {
#delimit ;                                      
coefplot 
(w21m`i', msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21c`i', msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w34m`i', msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21mf`i', msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cf`i', msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mf`i', msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mm`i', msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w21cm`i', msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w34mm`i', msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.005(0.0025)0.005, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( w21m w21c w34m  w21mf w21cf w34mf  w21mm w21cm w34mm )  coeflabels(w21m="M21" w21c="C21" w34m="M34"   w21mf="M21" w21cf="C21" w34mf="M34"  w21mm="M21" w21cm="C21" w34mm="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 8 "Men", color(black) size(large)) text(-0.0045 5 "Women", color(gray) size(large)) text(-0.0045 2 "All", color(red) size(large)) saving(`i', replace)   ;
#delimit cr
}


******************************************************************
*Figure A11: Alternative Specifications - Marital Status
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

gen w21mfM = workXtworkp_pc
gen w21cfM = workXtworkp_pc
gen w34mfM = workXtworkp_pc
gen w21mfS = workXtworkp_pc
gen w21cfS = workXtworkp_pc
gen w34mfS = workXtworkp_pc

quietly reg voted_munic  i.worker w21mfM   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1      & married==1   
quietly estimates store w21mfM
quietly reg voted_county i.worker w21cfM   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1      & married==1    
quietly estimates store w21cfM
quietly reg voted_munic  i.worker w34mfM   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1      & married==1    
quietly estimates store w34mfM
quietly reg voted_munic  i.worker w21mfS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1      & married==0     
quietly estimates store w21mfS
quietly reg voted_county i.worker w21cfS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1      & married==0    
quietly estimates store w21cfS
quietly reg voted_munic  i.worker w34mfS   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1      & married==0    
quietly estimates store w34mfS 

#delimit ;                                      
coefplot 
(w21mfM, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cfM, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mfM, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mfS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
(w21cfS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
(w34mfS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
, ylabel(-0.005(0.0025)0.0075, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)7)) 
keep(   w21mfM w21cfM w34mfM w21mfS w21cfS w34mfS    )  coeflabels(   w21mfM="M21" w21cfM="C21" w34mfM="M34"  w21mfS="M21" w21cfS="C21" w34mfS="M34"           )  
yline(0, lcolor(gray)) xline(3.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0045 2 "Married", color(gray) size(large)) text(-0.0045 5 "Single", color(gray) size(large)) ;
#delimit cr

 
******************************************************************
*Figure A12: Robustness to Dependents
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Dependents
#delimit ;
 twoway (scatter  tdepp_pc tworkp_pc if year==1921, mcolor(gray) m(oh)) 
 (scatter  tdepp_pc tworkp_pc if year==1934, mcolor(erose) m(oh)) 
 (lfit  tdepp_pc tworkp_pc if year==1921, lwidth(thick) lcolor(black))
 (lfit  tdepp_pc tworkp_pc if year==1934, lwidth(thick) lcolor(red)), 
 ytitle(Dependents Among Neighbors (%), size(vlarge)) ylabel(0(20)80, labsize(vlarge)) 
 xtitle(Workers Among Neighbors (%), size(vlarge)) xlabel(0(25)100, labsize(vlarge)) legend(off) scheme(s1color) ;
#delimit cr
 
*Sub-Figure a: Widows
#delimit ;
 twoway (scatter  twidp_pc tworkp_pc if year==1921, mcolor(gray) m(oh)) 
 (scatter  twidp_pc tworkp_pc if year==1934, mcolor(erose) m(oh)) 
 (lfit  twidp_pc tworkp_pc if year==1921, lwidth(thick) lcolor(black))
 (lfit  twidp_pc tworkp_pc if year==1934, lwidth(thick) lcolor(red)), 
 ytitle(Widows Among Neighbors (%), size(vlarge)) ylabel(0(20)80, labsize(vlarge)) 
 xtitle(Workers Among Neighbors (%), size(vlarge)) xlabel(0(25)100, labsize(vlarge)) legend(off) scheme(s1color) ;
#delimit cr 
 
 
******************************************************************
*Figure A13: Logit Models
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: Working Class
quietly logit voted_munic  i.worker w21m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
quietly estimates store w21m
quietly logit voted_county i.worker w21c   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
quietly estimates store w21c
quietly logit voted_munic  i.worker w34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
quietly estimates store w34m
quietly logit voted_munic  i.worker w21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
quietly estimates store w21mf
quietly logit voted_county i.worker w21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
quietly estimates store w21cf
quietly logit voted_munic  i.worker w34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
quietly estimates store w34mf 
quietly logit voted_munic  i.worker w21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
quietly estimates store w21mm
quietly logit voted_county i.worker w21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
quietly estimates store w21cm
quietly logit voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 
quietly estimates store w34mm
 
#delimit ;                                      
coefplot 
(w21m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21c, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(w21mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w21cm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(w34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.025(0.025)0.05, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( w21m w21c w34m  w21mf w21cf w34mf  w21mm w21cm w34mm )  coeflabels(w21m="M21" w21c="C21" w34m="M34"   w21mf="M21" w21cf="C21" w34mf="M34"  w21mm="M21" w21cm="C21" w34mm="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.023 8 "Men", color(black) size(large)) text(-0.023 5 "Women", color(gray) size(large)) text(-0.023 2 "All", color(red) size(large)) ;
#delimit cr

*Sub-Figure b: Upper and Middle Class
quietly logit voted_munic  i.upmid um21m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   
quietly estimates store um21m
quietly logit voted_county i.upmid um21c   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   
quietly estimates store um21c
quietly logit voted_munic  i.upmid um34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  
quietly estimates store um34m
quietly logit voted_munic  i.upmid um21mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1 
quietly estimates store um21mf
quietly logit voted_county i.upmid um21cf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1 
quietly estimates store um21cf
quietly logit voted_munic  i.upmid um34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1 
quietly estimates store um34mf 
quietly logit voted_munic  i.upmid um21mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==0  
quietly estimates store um21mm
quietly logit voted_county i.upmid um21cm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==0 
quietly estimates store um21cm
quietly logit voted_munic  i.upmid um34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0 
quietly estimates store um34mm
 
#delimit ;                                      
coefplot 
(um21m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um21c, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)   ) )
(um21mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um21cf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(um21mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(um21cm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
(um34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)   ) )
, ylabel(-0.025(0.025)0.05, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)10)) 
keep( um21m um21c um34m  um21mf um21cf um34mf  um21mm um21cm um34mm )  coeflabels(um21m="M21" um21c="C21" um34m="M34"   um21mf="M21" um21cf="C21" um34mf="M34"  um21mm="M21" um21cm="C21" um34mm="M34"  , angle(45)      )  
yline(0, lcolor(gray)) xline(3.5 6.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.023 8 "Men", color(black) size(large)) text(-0.023 5 "Women", color(gray) size(large)) text(-0.023 2 "All", color(red) size(large)) ;
#delimit cr
 
 
******************************************************************
*Figure A14: Alternative - New Voters
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*Sub-Figure a: working class
quietly reg voted_munic  i.worker w34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  & age<=27
quietly estimates store w34m
quietly reg voted_munic  i.worker w34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & age<=27  
quietly estimates store w34mf 
quietly reg voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & age<=27  
quietly estimates store w34mm
quietly reg voted_munic  i.worker w34m   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  & age>27
quietly estimates store w34mE
quietly reg voted_munic  i.worker w34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & age>27  
quietly estimates store w34mfE
quietly reg voted_munic  i.worker w34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & age>27  
quietly estimates store w34mmE

#delimit ;                                      
coefplot 
(w34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)) )
(w34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)) )
(w34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)) )
(w34mE, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick) lpattern(dash)) )
(w34mfE, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick) lpattern(dash)) )
(w34mmE, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick) lpattern(dash)) )
, ylabel(-0.005(0.005)0.015, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)4)) 
keep( w34m   w34mf   w34mm )  coeflabels( w34m="All"  w34mf="Women" w34mm="Men"  )   
yline(0, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) ;
#delimit cr

*Sub-Figure b: upper and middle class
quietly reg voted_munic  i.upmid um34m   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  & age<=27
quietly estimates store um34m
quietly reg voted_munic  i.upmid um34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & age<=27  
quietly estimates store um34mf 
quietly reg voted_munic  i.upmid um34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & age<=27  
quietly estimates store um34mm
quietly reg voted_munic  i.upmid um34m   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934  & age>27
quietly estimates store um34mE
quietly reg voted_munic  i.upmid um34mf   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1  & age>27  
quietly estimates store um34mfE
quietly reg voted_munic  i.upmid um34mm   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==0  & age>27  
quietly estimates store um34mmE

#delimit ;                                      
coefplot 
(um34m, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick)) )
(um34mf, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)) )
(um34mm, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick)) )
(um34mE, msym(o) mcolor(red) mlabcolor(red) msize(large) ciopts(lcolor(red) lwidth(midthick) lpattern(dash)) )
(um34mfE, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick) lpattern(dash)) )
(um34mmE, msym(o) mcolor(black) mlabcolor(black) msize(large) ciopts(lcolor(black) lwidth(midthick) lpattern(dash)) )
, ylabel(-0.005(0.005)0.015, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)4)) 
keep( um34m   um34mf   um34mm )  coeflabels( um34m="All"  um34mf="Women" um34mm="Men"   )  
yline(0, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) ;
#delimit cr
 
 
******************************************************************
*Figure A15: Alternative - Employment Heterogeneity
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

*sub-figure a: excluding domestic servants
gen w21mfD = workXtworkp_pc
gen w21cfD = workXtworkp_pc
gen w34mfD = workXtworkp_pc
gen w21mfDS = workXtworkp_pc
gen w21cfDS = workXtworkp_pc
gen w34mfDS = workXtworkp_pc

quietly reg voted_munic  i.worker w21mfD   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1       & dworker==0    
quietly estimates store w21mfD
quietly reg voted_county i.worker w21cfD   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1       & dworker==0    
quietly estimates store w21cfD
quietly reg voted_munic  i.worker w34mfD   C_age C_age2 i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1       & dworker==0    
quietly estimates store w34mfD 
quietly reg voted_munic  i.worker w21mfDS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1      & married==0 & dworker==0    
quietly estimates store w21mfDS
quietly reg voted_county i.worker w21cfDS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1      & married==0  & dworker==0    
quietly estimates store w21cfDS
quietly reg voted_munic  i.worker w34mfDS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1      & married==0 & dworker==0    
quietly estimates store w34mfDS 

#delimit ;                                      
coefplot 
(w21mfD, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cfD, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mfD, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mfDS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
(w21cfDS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
(w34mfDS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
, ylabel(-0.005(0.0025)0.0075, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)7)) 
keep(   w21mfD w21cfD w34mfD w21mfDS w21cfDS w34mfDS    )  coeflabels(   w21mfD="M21" w21cfD="C21" w34mfD="M34"  w21mfDS="M21" w21cfDS="C21" w34mfDS="M34"           )  
yline(0, lcolor(gray)) xline(3.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0038 2 "Excluding", color(gray) size(large)) text(-0.0038 5 "Excluding", color(gray) size(large)) 
text(-0.0045 2 "Domestic Workers", color(gray) size(large)) text(-0.0045 5.3 "Domestic W. & Married", color(gray) size(large)) ;
#delimit cr

*sub-figure b: excluding factory workers
gen w21mfF = workXtworkp_pc
gen w21cfF = workXtworkp_pc
gen w34mfF = workXtworkp_pc
gen w21mfFS = workXtworkp_pc
gen w21cfFS = workXtworkp_pc
gen w34mfFS = workXtworkp_pc

quietly reg voted_munic  i.worker w21mfF   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1       & fworker==0    
quietly estimates store w21mfF
quietly reg voted_county i.worker w21cfF   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1       & fworker==0    
quietly estimates store w21cfF
quietly reg voted_munic  i.worker w34mfF   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1       & fworker==0    
quietly estimates store w34mfF 
quietly reg voted_munic  i.worker w21mfFS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1921   & female==1      & married==0 & fworker==0    
quietly estimates store w21mfFS
quietly reg voted_county i.worker w21cfFS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_county==1 & year==1921   & female==1      & married==0  & fworker==0    
quietly estimates store w21cfFS
quietly reg voted_munic  i.worker w34mfFS   C_age C_age2  i.married   i.property_id, cluster(property_id), if right_munic==1  & year==1934   & female==1      & married==0 & fworker==0    
quietly estimates store w34mfFS 

#delimit ;                                      
coefplot 
(w21mfF, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21cfF, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w34mfF, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)   ) )
(w21mfFS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
(w21cfFS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
(w34mfFS, msym(o) mcolor(gray) mlabcolor(gray) msize(large) ciopts(lcolor(gray) lwidth(midthick)  lpattern(dash) ) )
, ylabel(-0.005(0.0025)0.0075, labsize(vlarge)) ytitle("Average Neighbor Effects", size(vlarge) )  xlabel(,labsize(vlarge) ) xscale(range(0(1)7)) 
keep(   w21mfF w21cfF w34mfF w21mfFS w21cfFS w34mfFS    )  coeflabels(   w21mfF="M21" w21cfF="C21" w34mfF="M34"  w21mfFS="M21" w21cfFS="C21" w34mfFS="M34"           )  
yline(0, lcolor(gray)) xline(3.5, lcolor(gray)) levels(95) scheme(s1mono) legend(off) vertical title("", size(vlarge) ) ysize(5) 
text(-0.0038 2 "Excluding", color(gray) size(large)) text(-0.0038 5 "Excluding", color(gray) size(large)) 
text(-0.0045 2 "Factory Workers", color(gray) size(large)) text(-0.0045 5.3 "Factory W. & Married", color(gray) size(large)) ;
#delimit cr


******************************************************************
*Table A11: Alternative - Proxy Voting
******************************************************************
use dta\sod2134, clear                       //using sod2134 data set

reg voted_munic   i.class_wife C_age C_age2  i.married##i.female   i.property_id ,cluster(property_id) , if year==1921 & right_munic==1 
reg voted_county  i.class_wife C_age C_age2  i.married##i.female   i.property_id ,cluster(property_id) , if year==1921 & right_county==1 
reg voted_munic   i.class_wife C_age C_age2  i.married##i.female   i.property_id ,cluster(property_id) , if year==1934 & right_munic==1 
 
reg voted_munic    C_age C_age2  i.married##i.female   i.property_id ,cluster(property_id) , if year==1921 & right_munic==1  & worker==1
reg voted_county   C_age C_age2  i.married##i.female   i.property_id ,cluster(property_id) , if year==1921 & right_county==1 & worker==1
reg voted_munic    C_age C_age2  i.married##i.female   i.property_id ,cluster(property_id) , if year==1934 & right_munic==1  & worker==1
  
  
******************************************************************
*Figure A16: Sweden in Context
******************************************************************

*sub-figure a: female labour force
use dta\lab, clear                       //using labour data set
#delimit ;
twoway (scatter female_labour lab_lbf, mcolor(black) m(O) mlab(female_labour) mlabpos(12)) 
, ytitle(Female Labour Force (%), size(med)) ylabel(0(25)75, labsize(med)) 
xtitle("", size(vlarge)) xlabel(1(1)13, valuelabel angle(45) labsize(med))
title("", nobox) note("", nobox) legend(off) scheme(s1mono) ysize(4) ;
#delimit cr

*sub-figure b: industrial workers
use dta\lab, clear                       //using labour data set
#delimit ;
twoway (scatter ind_workers lab_ind, mcolor(black) m(O) mlab(ind_workers) mlabpos(12)) 
, ytitle(Industry Workers (%), size(med)) ylabel(0(25)75, labsize(med)) 
xtitle("", size(vlarge)) xlabel(1(1)13, valuelabel angle(45) labsize(med))
title("", nobox) note("", nobox) legend(off) scheme(s1mono) ysize(4) ;
#delimit cr















